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Abstract. We present BeppoSAX results of a spatially resolved spectral analysis of A3571, a relaxed nearby cluster 
of galaxies. In the central 2' (130 h^ kpc) radius the metal abundance is 0.49 ± 0.08 solar and the absorption 
(1.13±0.28) x 10 21 atom cm 2 , whereas elsewhere within an 8' (520 hgQ 1 kpc) radius the abundance is 0.32 ±0.05 
solar and the absorption consistent with the galactic value of 4.4 x 10 20 atom cm -2 . The significant central metal 
abundance enhancement is consistent with the supernova enrichment scenario. The excess absorption may be 
attributed to the cooling flow, whose mass flow rate is 80 ± 40 Mq yr from our spectral fit. The BeppoSAX 
and ASCA radial temperature profiles agree over the entire overlapping radial range r < 25' = 1.6 h^" Mpc. The 
combined BeppoSAX and ASCA temperature profile exhibits a constant value out to a radius of ~10' (650 h^" kpc) 
and a significant decrease (T oc r _0 ' 55 i corresponding to 7 = 1.28) at larger radii. These temperature data are 
used to derive the total mass profile. The best fit NFW dark matter density model results in a temperature profile 
that is not convectively stable, but the model is acceptable within the uncertainties of the data. The temperature 
profile is acceptably modeled with a "core" model for the dark matter density, consisting of a core radius with a 
constant slope at larger radii. With this model the total mass and formal 90% confidence errors within the virial 
radius ri78 (2.5 h^, 1 Mpc) are 9.ll 3 'g x 10 14 h^, 1 Mq, by a factor of 1.4 smaller than the isothermal value. The 

, gas mass fraction increases with radius, reaching fgas(ri7s) = 0.26^q'5o X hj . Assuming that the measured gas 

mass fraction is the lower limit to the primordial baryonic fraction gives f2 m < 0.4 at 90% confidence. 
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1. Introduction 

Abell 3571 (hereafter A3571) is a nearby (z=0.04) bright 
and hot (~ 7 keV) cluster of galaxies, situated in the 
Shapley supercluster. A3571 is a rich (> 100 galaxies) 
Bautz- Morgan type I cluster (Abell et al. 1989). It is 
known to possess at its center a giant galaxy MCG05-33- 
002 with extensive optical halo (Kemp & Meaburn 1991). 
The galaxy distribution, as well as the optical halo, are 
aligned in the direction of the major axis of the central 
galaxy. Even though the optical data imply dynamical ac- 
tivity, the X-ray observations have revealed that A3571 is 
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a relaxed cluster. The ROSAT PSPC (position sensitive 
proportional counter) image of this cluster shows no sub- 
structure and no deviations from the azimuthal symmetry 
(Nevalainen et al. 2000) while the ASCA temperature map 
shows no azimuthal structure indicative of cluster mergers 
(Markevitch et al. 1998). A3571 has a considerable cool- 
ing flow (Peres et al. 1998), also indicative of a relaxed 
cluster. 

In this paper we report on a BeppoSAX observation 
of A3571. We use the data to perform spatially resolved 
spectroscopy in order to determine the temperature pro- 
file of A3571. Using the BeppoSAX temperature profile 
derived here, together with the ROSAT brightness profile 
(Nevalainen et al. 2000), and ASCA temperature profile 
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of A3571 (Markevitch et al. 1998), we constrain the total 
mass in A3571, assuming hydrostatic equilibrium. 

We use H = 50 h 50 km s" 1 Mpc~\ m = 1, A = 0, 
which corresponds to the linear scale of 65 kpc/arcmin at 
the cluster redshift z = 0.04. We consider uncertainties 
and significances at 90% level throughout the paper. 

2. Observations 

The cluster of galaxies A3571 was observed by BeppoSAX 
(Boella et al. 1997a) between 2000 February 2 17:42 and 
February 6 07:41 UTC. In this paper we discuss data from 
the imaging Medium-Energy Concentrator Spectrometer 
(MECS; 1.8-10 keV; Boella et al. 1997b), and the imag- 
ing Low-Energy Concentrator Spectrometer (LECS; 0.1- 
10 keV; Parmar et al. 1997). The MECS consists of two 
identical grazing incidence telescopes with imaging gas 
scintillation proportional counters in their focal planes. 
The LECS uses an identical concentrator system to the 
MECS, but utilizes an ultra-thin entrance window and a 
driftless configuration to extend the low-energy response 
to 0.1 keV. The fields of view (FOV) of the LECS and 
MECS are circular with diameters of 37' and 56', respec- 
tively. The MECS has a 2' thick circular window support 
structure, or strongback, which is centered 10' from the 
FOV center. In addition, four 2' thick radial spokes extend 
outwards from the circular support. The LECS strongback 
consists of a square grid with a separation of 4'. Standard 
screening and reduction procedures were applied to pro- 
duce linearized event files using the SAXDAS 2.0.0 pack- 
age. The exposure times for the MECS and LECS were 
66 ks and 28 ks, respectively. 

3. Spatially resolved spectroscopy 

The energy dependent Point Spread Functions (PSF) of 
the LECS and MECS instruments complicate the analy- 
sis of extended objects such as clusters of galaxies. In the 
case of the MECS the problem is less severe, because the 
mirror and the detector PSF partially cancel each other. 
For the LECS, the 90% encircled energy radius increases 
from 6' to 9' between 8 keV and 0.2 keV. In this work, 
we accounted for the energy dependent PSF by generat- 
ing appropriate instrument response files to be used when 
fitting the accumulated spectra. For the response gener- 
ation, we used the ROSAT brightness profile of A3571 
(Nevalainen et al. 2000) to model the flux distribution. 
For the MECS we produced corrected auxiliary responses 
with the SAXDAS Effarea tool, as was used for the clus- 
ters A2319 (Molendi et al. 1999) and A2256 (Molendi et 
al. 2000). Effarea convolves the brightness model with an 
analytical model for the MECS PSF to calculate the in- 
strumental response (see D'Acri et al. 1998 for details). 
Effarea corrects for the energy-dependent vignetting. For 



the LECS we used the ray-tracing code LEMAT (Lammers 
1997) to produce appropriate responses, taking into ac- 
count the extraction region size and source position on 
the detector, mirror vignetting and obscuration by the 
strongback. For any given extraction region and energy 
LEMAT estimates the contribution of photons from re- 
gions external to the extraction region, and the number of 
photons originating from within the extraction region that 
are detected externally (see Kaastra et al. (1999) for the 
application of LEMAT in the case of the A2199 cluster). 

We accumulated LECS and MECS spectra in annular 
regions centered on the X-ray emission maximum in the 
MECS (see Fig. |l|). The bounding radii for MECS and 
LECS spectra are 0'-2', 2'-4', 4'-6', 6'-8', and 12'-16'. 
Between 16'-25' only MECS data were used due to the 
smaller LECS FOV. The MECS data between radii of 8'- 
12' are shadowed by the strongback, and the LECS data 
alone do not give adequate constraints on the spectral pa- 
rameters. Between radii of 12'-16' the regions shadowed 
by the strongback are excluded from the MECS spectra. 
Between 16'-25' the regions affected by the strongback 
and the calibration sources are excluded. We obtained 
background spectra from standard blank sky observations 
using the same extraction regions as for the source spec- 
tra. LECS and MECS data were selected in the energy 
ranges 0.1-4.0 keV and 1.8-10 keV, respectively. 

The LECS and MECS spectra were rebinned to over- 
sample the full-width half maximum of the energy reso- 
lution by a factor 3 and additionally to have a minimum 
of 20 counts per bin to allow use of the \ 2 statistic, ex- 
cept at radii of 12'-16' and 16'-25' where a minimum of 
30 and 40 counts per bin was used to improve the signal 
to noise ratio. At given radii (below 16') the MECS and 
LECS spectra were fit simultaneously allowing the rela- 
tive normalization between the two instruments to vary 
by ±30% to allow for cross-calibration uncertainties (e.g., 
Fiore et al. 1999). We allowed a 2% systematic error in 
the fits to account for absolute calibration uncertainties. 

3.1. Radial analysis 
3.1.1. Global fit 

To obtain the global properties of the intracluster gas, we 
fitted the MECS and LECS data between radii 2'-8'. The 
inner 2' region was excluded from the fit, because A3571 
is known to have a moderate cooling flow in the center 
(see below) , and the data beyond 8' were excluded due to 
strongback obscuration. We model the emission from the 
hot intracluster gas using the MEKAL model (Mewe et al. 
1995) in XSPEC. A single temperature fit to 2-8' spec- 
tra with the absorption, Nh fixed at the galactic value of 
4.4 x 10 20 atom cm~ 2 (Murphy et al. 2000) is unaccept- 
able with a x 2 of 120.0 for 95 degrees of freedom (dof). 
There are systematic residuals below 0.3 keV (positive) 
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and between 0.3-0.6 keV (negative), (see Fig. ||). These 
features suggest the presence of a warm absorber, which 
has been detected in several cooling flow clusters (Buote 
2000a, 2000b). To get an order-of-magnitude estimate of 
the properties of this component, we added an absorption 
edge to the model to account for the extra absorption be- 
tween 0.3-0.5 keV and an additional mekal component 
to account for the excess emission below 0.3 keV. We kept 
Nh fixed at the galactic value and constrained the metal 
abundances of both thermal components to be equal, since 
the data are not adequate to determine them both sepa- 
rately. The resulting fit is acceptable with a % 2 /dof of 
90.3/91 and requires an absorption depth of 5. 9~_2 7 & t an 
edge energy of 0.27 ± 0.01 keV. The warm component has 
T = 0.18 ± 0.03 keV and a normalization 0.27t£; 2 ^ that 
of the hot component. 

However, this model is not physically meaningful, be- 
cause the edge energy and the temperature are not linked 
by physical processes. Therefore we built a warm ab- 
sorber model by computing the transmission of the MEKAL 
model as a function of energy for a grid of temperatures, 
abundances and column densities, and imported this table 
model into XSPEC. We fitted the data with a model con- 
sisting of the above warm absorber model, Galactic Nh, 
a mekal component and an additional MEKAL compo- 
nent, whose temperature we fixed to that of the absorber. 
We were not able to explain the observed features with 
this model. To produce the observed depletion at 0.3-0.5 
keV the absorber should have a temperature below ^0.01 
keV, but at such low temperatures the model does not 
produce the observed excess emission at 0.1-0.3 keV. Also, 
the oxygen absorption edge at 0.6 keV, predicted with this 
model, is not observed. The observed features may be due 
to LECS background or response uncertainties, or due to a 
more complex absorber. At the lowest energies the back- 
ground level is comparable to the cluster emission. The 
PSF correction applied to the LECS response is largest 
at the lowest energies and so the uncertainties may be 
stronger there. Also, the galactic absorption reduces the 
intrinsic flux by three orders of magnitude at 0.1 keV, and 
small deviations in the abundances of the absorbing ma- 
terial could produce low-energy structure. More complex 
absorption could result from ionized or partially covered 
material. However, the quality of the data is insufficient 
to disentangle these effects. 

Fitting the data with a model consisting of a Galactic 
Nh and a single mekal component, excluding bins below 
0.6 keV leads to an acceptable fit (see Table 1). As seen in 
Fig. ||, the slight miscalibration of the MECS pulseheight- 
channel to energy relation results in slight offset between 
the iron 6.7 keV line of the plasma model and the data (as 
noticed in the case of A2256 by Molendi et al. 2000). The 
reason for this is unclear. The in-flight response calibra- 
tion is made using Crab data, and possible gain problems 



should have been taken care of by the response creation 
procedures. At the moment it is not possible to accurately 
correct for this effect. However, the mismatch of the iron 
line energy has a negligible effect on temperatures. We 
verified this by fitting the data excluding the 6-7 keV 
band and noticing that the changes in best-fit tempera- 
tures were less than the parameter uncertainties. We esti- 
mated the effect of the gain miscalibration on the tempera- 
ture measurements by adjusting the gain linear coefficient 
by 1%, after which the model iron line matched the data, 
and redoing the fits with this adjustment (see Table 1). 
This does not properly correct for the gain problem, but 
we use this as an estimate of the systematic uncertainty 
due to gain miscalibration. This adjustement results in 
systematic shift towards lower temperatures of about the 
size of the parameter uncertainties at radii below 8' . At 
larger radii the effect is negligible. In the following, when 
presenting results we include this uncertainty in the pa- 
rameter values. 

3.1.2. Center 

In the 0'-2' radial bin the spectral fit is formally accept- 
able with a model consisting of a single absorber and a 
single temperature plasma component with Galactic Nh- 
However, the fit improves significantly, when Nh is al- 
lowed to be a free parameter. The resulting absorption 
((11.3^2 6) x 10 2 ° atom cm -2 ) is significantly (by a fac- 
tor of 2.6) above the galactic value. Qualitatively simi- 
lar behavior has been found in other cooling flow clusters 
(e.g. Allen & Fabian 1997; Molendi et al. 1998; Molendi 
& De Grandi, 1999; De Grandi & Molendi 1999b). The 
corresponding BeppoSAX metal abundance (0.49±0.08 in 
solar units) is significantly (by a factor of 1.5) above the 
global value. A study of several clusters, including A3571 
using ASCA data (Dupke & White 2000) also shows cen- 
tral metal enhancement and gives for A3571 abundance 
values consistent with ours. 

The cooling flow found in the center of A3571 (Peres 
et al. 1998) contributes a comparable amount of flux to 
the ambient gas in the central bin. To study this more we 
added a cooling flow model (mkcflow in XSPEC) fixing 
the lower boundary of the continous temperature range to 

0. 1 keV and forcing the upper boundary to equal the am- 
bient temperature. The metal abundance was constrained 
to be equal to the ambient one, because the data qual- 
ity does not allow them to be measured separately. We 
also fixed the mass flow rate to 75 M Q yr" 1 as found in 
ROSAT analysis (Peres et al. 1998). As a result, the tem- 
perature of the ambient gas in the central bin becomes 
9.5^'g keV. Letting the mass flow rate to be a free param- 
eter, the temperature of the gas becomes very uncertain, 
but nevertheless the cooling flow is significantly detected, 

1. e. M = 80 ± 40 Moyr" 1 at 90% confidence level. 
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Table 1. Fit results. The model used in the fits is mekal absorbed by Galactic Nh, unless stated differently in remarks 
column. The errors are given at 90% confidence level 



NO GAIN CORRECTION GAIN CORRECTION 
Radius T Abund xVdof T Abund X 2 /dof Remarks 
(keV) (solar) (keV) 



2-8 


7 q+o.4 


n OO+0.05 
u ' oz -0.04 


98.5/82 


7 6 +0 - 3 

'• D -0.3 


31 + 05 

u - ol -0.03 


95.2/81 




0-2 


o 1+0.5 
8-1-0.4 


o ™ +0 - 08 

u - ou -0.07 


87.7/82 


7 y+0.5 

'•'-0.3 


n 4Q+ 009 

u -* y -0.07 


87.7/82 


wabs x mekal, Nh = Galactic 


0-2 


'• D -0.5 


4Q+ 08 
u '^ y -0.07 


66.8/81 


7 O+0.4 
'• Z -0.5 




61.4/81 


wabs x mekal, Nh free 


0-2 


9 5+ 11 


n co+0.09 
U ' JO -0.08 


56.1/81 


o y+0.9 

°-'-0.7 


u - dz -0.08 


53.1/81 


wabs x (mekal + mkcffow), Nh free 


2-4 


7 6 +0 ' 5 

'• D -0.3 


31 + 06 


97.8/83 


7 O+0.4 

'•■J-o.s 


n on+0.05 
u -° -0.06 


92.2/83 




4-6 


7.8±°i 


S7+ 07 

u -°' -0.07 


95.3/81 


74 +o.5 

'•^-0.4 


37+ 07 

u -°' -0.07 


92.9/81 




6-8 


8.41°;? 


26+ 09 


66.4/82 


8 .oir 6 


0-25lH 9 9 


68.7/82 




12-16 


6-6±i; 


u - oz -0.20 


75.8/71 


6-4±i; 


0.29±g;?S 


76.3/71 




16-25 




n is +0 ' 37 

u - 1,5 -0.13 


61.2/47 


4.2±J:S 


o.i9±S:S 


60.8/47 





3.1.3. Other radii 

Using a single mekal model (as explained above), we find 
that temperatures and metal abundances in the radial bins 
2'-4', 4'-6' and 6'-8' are consistent with being constant 
(see Fig. | for the T profile). Between 12'-16' the best 
fit temperature is lower, but still agrees within the uncer- 
tainties with the 2'-8' fit results. In the 16'-25' bin the 
temperature is significantly lower than in the other bins. 
The BeppoSAX temperature profile of A3571 derived here 
has slightly higher values than that derived using ASCA 
for this cluster (Markevitch et al. 1998), but both agree 
within 90% confidence uncertainties at each radius. 

To quantify the behavior of the temperature profile, 
we modeled the 3-dimensional temperature profile with a 
constant value T const out to r = 10', and with a behav- 
ior T oc r _1 at larger radii. We projected this model into 
BeppoSAX and ASCA bins and iteratively determined 
the best-fit parameters to the data and to the Monte - 
Carlo sets (see below) in order to obtain the uncertainties. 
We obtain T const = 7.8±o'J and i= 0.55±q%1 which corre- 
sponds to a polytropic model T oc pj^, 1 with 7 = 1.281q^ 2 
beyond r = 10'. However, because the observed T is con- 
stant up to r = 10' and therefore different from the gas 
density profile, the polytropic model does not give an ac- 
ceptable fit to the data. 

Due to the strongback obscuration and low count rates 
at large radii the BeppoSAX data are not ideal for mea- 
suring the spatial variation of the spectral parameters in 
A3571. However, we divided the 2'-8' MECS and LECS 
data into four sectors and between radii 12' and 25' we 
used two sectors of MECS data. The temperatures and 
metal abundances in these sectors are consistent with the 
average values, i.e. we find no spatial variation. This is 
consistent with the ASCA results of A3571 (Markevitch 
et al. 1998). 



3.1.4. Fluxes 

With the best-fit models (using the MECS normal- 
izations) we determined the fluxes and luminosities 
in complete annuli in different radial bins. Above 8', 
where the sky coverage is incomplete due to exclusion 
of strongback and calibration source regions, we assume 
azimuthal symmetry when computing the fluxes. In the 
2'-8' region we measure unabsorbed 2-10 keV fluxes and 
luminosities of F2-10 = 7.3 ± 0.4 10 -11 erg s -1 cm" 2 
and L2-10 = 5.1 ± 0.2 hj~ 2 10 44 erg s" 1 and bolomet- 
ric values of Fboi = 15.5 ± 0.7 10" 11 erg s" 1 cm" 2 , 
L bo i = 10.8 ± 0.4 h^ 2 10 44 erg s" 1 . For compari- 
son, we distribute the measured ROSAT flux of 
A3571 (Markevitch 1998) between different radii 
using the best-fit (3 profile from Nevalainen et al. 
(2000). Between 2'-8' the 0.1-2.4 keV ROSAT flux 
is 5.9 ± 0.3 10" 11 erg s" 1 cm~ 2 and agrees with the 
BeppoSAX value of 5.8 ± 0.1 10~ n erg s" 1 cm" 2 . 
Also within the 12'-16' annulus the ROSAT value of 
0.871q;q 4 10~ n erg s~ 4 cm" 2 is consistent with our value 
of 0.84±g;^ 10" 11 erg s" 1 cm" 2 . At the largest radii 
there is some model dependent discrepancy, because 
the ROSAT flux estimate is obtained with a single 
temperature model, whereas the BeppoSAX temperature 
decreases at large radii. However, the ROSAT value for 
radii of 16'-25' (0.96t°;°| 10" 11 erg s" 1 cm" 2 ) agrees 
with our value (l.Ctgj 10 -11 erg s _1 cm" 2 ) within 3cr 
confidence. 

4. Mass determination 

4.1. Method 

The ASCA temperature map (Markevitch et al. 1998) 
and the ROSAT brightness distribution (Nevalainen et al. 
2000) of A3571 show no deviations from azimuthal sym- 
metry and no substructure, and the BeppoSAX data do 
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not contradict this view. Therefore we are justified in as- 
suming that the hot X-ray emitting material in A3571 is 
in hydrostatic equilibrium. 

Since the BeppoSAX temperature profile agrees with 
the ASCA profile, we use both simultaneously to con- 
strain the total mass. For the details of the method we 
refer to the similar analysis of the ASCA data of A3571 
(Nevalainen et al. 2000) and A401 (Nevalainen et al. 1999). 
Briefly, we model the dark matter density with two differ- 
ent types of profiles. In the "core" model 



Pdark OC 1 



-a/2 



(1) 



the density is relatively constant inside the core radius and 
has a constant slope a at large radii, analoguous to the (3 
model for the gas density. The "cusp" profile 



r \ v / r 

Pdark OC I H 

\a d J V a d 



(2) 



has a changing slope and a central cusp. With r\ = 1 and 
a = 3, the cusp model corresponds to the "universal den- 
sity profile" which is a good description of the cluster cold 
dark matter halos in simulations of hierarchical cluster- 
ing (Navarro et al. 1995, 1997, hereafter NFW). Our data 
are not adequate to determine the shape of the cusp, and 
therefore we fix the shape parameter 77 to be 1 in the cusp 
models, as suggested by the above simulations, but vary 
the other parameters. We solve the hydrostatic equilib- 
rium equation (e.g., Sarazin 1988): 



M tot (< r) = 3.70 x 10 13 M ( 
dlnp ga 



T(r) r 
keV Mp^ 



d In ) 



rflnT 
rflnr 



(3) 



(using fi = 0.60), with respect to T, in terms of the dark 
matter and gas density profile parameters (the analytical 
solutions are given in Markevitch & Vikhlinin 1997). We 
fix the gas density to that found from the ROSAT anal- 
ysis (i.e. (3 model with slope parameter (3 = 0.68 ± 0.03 
and core radius & x = 3.85 ± 0.35' which describes well the 
PSPC surface brightness profile data in the radial range 
considered in this work (Nevalainen et al. 2000), calculate 
the 3-dimensional temperature profile model correspond- 
ing to given dark matter parameters, project it on the 
BeppoSAX and ASCA annuli, compare these values to 
the observed temperatures and iteratively determine the 
dark matter distribution parameters. 

To propagate the errors of the temperature profile data 
to our mass values, we repeat the procedure for a large 
number of Monte - Carlo temperature profiles with added 
random errors. The gain calibration uncertainty will only 
affect the BeppoSAX temperatures and its effect on the 
combined BeppoSAX + ASCA mass profile is negligible 



and we thus ignore it. From the distribution of the ac- 
ceptable Monte - Carlo models, we determine the la con- 
fidence intervals of the mass values as a function of radius. 
We convert these values to 90% confidence values, assum- 
ing a Gaussian probability distribution. 

4.2. Results 

Even though our BeppoSAX temperature profile is consis- 
tent with the ASCA profile (Markevitch et al. 1998), the 
BeppoSAX profile is somewhat steeper. As a consequence 
the best-fit cusp model becomes too steep to be convec- 
tively stable above r = 26', i.e. the effective polytropic 
index 



7 (r) 



d log T(r) 
d log p(r) 



1 



(4) 



exceeds §. However, some of the Monte 



Carlo models, 

which quantify the effect of the model parameter uncer- 
tainties, are convectively stable. Therefore the cusp model, 
as well as the NFW model are only marginally acceptable, 
and thus we report no mass values based on this model. 

The best-fit core model temperature profile is accept- 
able with a x 2 of 8.5 for 4 dof. This model is on the border 
of being convectively unstable, but can still be accepted 
(see Fig. || for the profile). The best-fit is obtained with 
a core radius of ad = lO.'l, a slope of a = 6.0, a cen- 
tral temperature of To = 6.2 keV and central dark matter 
density of pd(0) = 6.6 x 10 -26 g cm -3 (6 times as high 
as the central gas density). The temperature data covers 
the radial range up to 35'. Extrapolating the temperature 
profile slightly (to 39') using our best fit model, we can 
determine the cluster mass at the cosmologically interest- 
ing radius, rn$, i.e. the virial radius in (fi m — 1, J7a = 0) 
universe based on the spherical collapse model (Eke et al. 
1996). We obtain 



ri78 



38.4' = 2.5 h-J Mpc 



and the total mass within this radius is 



M 



ITS 



9.lt?;jjx 1Q 14 h$Q Mq, 



(5) 



(6) 



(see Table g for the mass values at several interesting 
radii). We note that half of the Monte - Carlo fits used for 
the mass uncertainty estimation are convectively unsta- 
ble. Models with larger negative temperature gradients at 
a given radius are more massive if the temperatures at that 
radius are equal (see Eq. ||). Thus, applying the convec- 
tive stability constraint would produce a strong constraint 
for the upper bound to the mass profile (e.g. at ins we 
would obtain Mi 78 = 9.1±?"| x 10 14 M Q ). However, 
such a constraint is model dependent and does not allow 
the uncertainties in the temperature measurements to be 
propagated to the mass estimates. We therefore choose not 
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Table 2. Mass determinations for A3571 assuming 
H = 50 h 50 km s _1 Mpc -1 . All uncertainties are formal 
90% confidence limits resulting from our analysis 



Parameter 



Value 



M tot (< 3.85' = a x ) [M ] 
M tot (< 15.4' = 1 Mpc ) [Mq] 
M to t(< 25.9' =r 50 o) [M ] 
M tot (< 35' ) [M ] 
M tot (< 38.4' =ri 78 ) [M ] 

f gas (< 25.9' =r 50 o) x h.% 2 



5.9 



hl.7 
-1.2 
K Q+0.8 
°- y -0.8 

7.9±H x 10 



x 10 13 
x 10 14 

,14 



8.8: 



1.5 



x 10 1 



9.i x io 14 



fgas(< 38.4' 



n , 3/2 

ri7s) x h 5 ^ 



0.19 
0.26 



+0.03 
0.06 

+0.05 
0.10 



to use this constraint to reject any of the Monte - Carlo 
fits. 

The mass values obtained here are consistent with 
the corresponding values obtained using only ASCA 
(Nevalainen et al. 2000), and more precise. In addition to 
the better data, the improvement in the precision is due to 
the fact that in the present work we cannot apply the cusp 
model (due to convective instability) as in Nevalainen et 
al. (2000). However, if we were to apply the convective sta- 
bility constraint for the Monte - Carlo fits with the core 
model, the mass uncertainties derived here would be much 
smaller. Our mass values within rsoo and rns are smaller 
by factors of 1.1 and 1.4, respectively, than the mass ob- 
tained assuming that the cluster is isothermal with the 
emission-weighted ASCA value Tx — 6.9 keV (Markevitch 
1998) but consistent within the uncertainties. Within r 50 o 
our total mass value is significantly smaller (by a factor 
of 1.6) than the value predicted by the cosmological sim- 
ulations of Evrard et al. (1996) for a cluster with T = 6.9 
keV. 

Our results indicate that the dark matter density falls 
as r~ 6 whereas the gas density falls as r~ 2 and thus the 
gas mass fraction increases rapidly at large radii. At rsoo 
our gas mass fraction value (see Table 2) is consistent 
with those for A2256 (Markevitch & Vikhlinin 1997), A401 
(Nevalainen et al. 1999), A2199 and A496 (Markevitch et 
al. 1999) and for samples in Ettori & Fabian (1999) and 
Mohr et al. (1999). At the virial radius we obtain 



f gas {< rns) 



2fi +0 05 h~ 

u - zo -0.10 rl 5 



3/2 



(7) 



Assuming that this is a lower limit on the cosmological 
baryon fraction, or that 



£2m 



Mb 

Mtnt 



> fgas(< n 78 ), 



(8) 



and using Q h = 0.076 ± 0.007 (Buries & Tytler 1998) and 
Ho > 50 km s^ 1 Mpc -1 we obtain an upper limit for the 
cosmological mass density parameter of f2 m < 0.41. 



5. Conclusions and discussion 

The heavy element abundace in the central 2' radial bin 
is significantly above the cluster average. Our values are 
consistent with the values obtained from ASCA data of 
A3571 (Dupke & White III 2000a) who also found sim- 
ilar enhancements in clusters A2199 and A496 (Dupke 
& White III 2000b). Also BcppoSAX analyses on several 
clusters (Irwin & Bregman 2000; Molendi & De Grandi 
1999; Dc Grandi & Molendi 1999b; Molendi et al. 2000) 
have yielded central metal enhancements. These obser- 
vations are consistent with the hypothesis of the metal 
enrichment of the intracluster gas by supernovae (e.g. 
Arimoto & Yoshii 1987; Matteucci & Gibson 1995). The 
central BeppoSAX absorption is significantly above the 
galactic value. Qualitatively similar behavior has been 
found in other cooling flow clusters (e.g. Allen & Fabian 
1997; Molendi et al. 1998; Molendi & De Grandi, 1999; De 
Grandi & Molendi 1999b). This suggests that the excess 
absorption is due to the cooling material accumulated in 
the central region of the cluster. 

The significantly decreasing temperature profile of 
A3571 derived here with BeppoSAX data agrees with the 
ASCA temperature profile of Markevitch et al. (1998). 
Similarly, the BeppoSAX temperature profiles for clusters 
A2029 (Molendi & De Grandi, 1999), A2256 (Molendi et 
al. 2000), and A3266 (De Grandi & Molendi 1999) show 
a temperature decline in good agreement with the respec- 
tive ASCA results (Sarazin et al., 1999; Markevitch & 
Vikhlinin, 1997; Markevitch et al., 1998), although the 
BeppoSAX profile for A2319 (Molendi et al. 1999) does 
not indicate a decline reported by Markevitch (1996) from 
ASCA data, even though both give consistent tempera- 
tures in the entire overlapping radial range. We used both 
BcppoSAX and ASCA temperature profiles to constrain 
the total mass of A3571, assuming hydrostatic equilib- 
rium. The BcppoSAX profile, even though it is formally 
consistent with the ASCA one, is somewhat steeper and 
results in problems with convective stability. The best fit 
cusp model (and the NFW universal profile) for the dark 
matter density are too steep compared to the gas density 
profile for the system to be convectively stable, but they 
are acceptable within the uncertainties of the data. The 
almost constant temperature up to r = 10' and steep de- 
crease at larger radii, as exhibited with both BeppoSAX 
and ASCA analysis cannot be acceptably modeled with 
polytropic models. The data are well modeled using the 
"core" model for the dark matter density, which is analo- 
gous to the (3 model for the gas density. The best-fit core 
model is at the border of being convectively unstable, but 
can still be accepted. The measured masses within rsoo 
and ri7§ are factors of 1.1 and 1.4 smaller than the val- 
ues obtained with the isothermal analysis of A3571 data, 
respectively, where the difference is consistent with the 
uncertainties. Within r 50 o the measured mass is signifi- 
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cantly smaller, by a factor of 1.6, than that predicted by 
the simulations of Evrard et al. (1996). 

The dark matter density profile in A3571 de- 
creases as r~ 6 , much faster than the gas density. 
Hence the gas mass fraction increases with radius, with 



f gas (ri78) = 0.26; 



-0.05 x h 3/2 



0.10 



Mo 



Assuming that this is the 



lower limit for the primordial baryonic fraction, we ob- 
tain Q m < 0.4 at 90% confidence. 
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RA [decimol degrees (J2000)] 

Fig. 1. BcppoSAX MECS 1.8-10 keV surface brightness image of A3571. The thick lines are the brightness contours 
after smoothing with a gaussian of a = V at values of 83, 33, 8, 3, 1 and 0.8% of the peak value of 1.2xl0~ 2 c s^ 1 
arcmin~ 2 . The thin lines show the regions where spectra were accumulated. The shaded area was excluded from the 
accumulations due to obscuration from the strongback 
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Fig. 2. The LECS and MECS spectra accumulated between 2'-8' together with the best-fit model single component 
(mekal) model absorbed by the galactic column. Significant structure is visible at energies <0.6 keV 
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Fig. 3. The solid crosses show the best-fit A3571 temperature with 90% confidence uncertainties obtained from 
BeppoSAX data using the MEKAL model. The error bars are expanded to include the shift due to the gain cali- 
bration uncertainty (i.e. the uncertainties include the 90% confidence intervals of the values in columns 2 and 5 at 
Table 1). At radii >2' the absorption consists of a galactic absorber. In the fit to the central 2' data, Nh was allowed 
to vary and a cooling flow component was included. The dotted line shows the single temperature fit. The dashed 
crosses show the values obtained using ASCA by Markevitch et al. (1998) 
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Fig. 4. The best-fit values and 90% confidence uncertainties of the metal abundance corresponding to temperature 
values in Fig. ||. The dashed lines show the global values obtained from the fit to the 2'-8' spectra 
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Fig. 5. The solid lines in panels a and b show the 3-dimcnsional best-fit core model temperature. The thick dotted 
lines in panels a and b show the best-fit core model temperatures projected onto the BeppoSAX and ASCA radial 
bins. The crosses are the temperature data as in Fig. 3. The thin dotted lines show a representative sample of Monte 
- Carlo models. Panel c shows the corresponding dark matter density profiles (solid and dotted lines) and the gas 
density profile (dashed line) 
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Fig. 6. The left panel shows the total mass profile (solid line) with 90% confidence uncertainties (shaded area) obtained 
with the core model. The dotted line shows the total mass obtained assuming a constant temperature T = 6.9 keV 
along the cluster. The dashed line shows the gas mass profile. The right panel shows the corresponding gas mass 
fraction profiles. Masses are evaluated using H = 50 km s" 1 Mpc -1 (total mass scales as H -1 and gas mass as H~ 5 / 2 ) 



